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Time-dependent  wave  packet  quantum  scattering  (TWQS)  calculations  are  presented  for 
HD^(i;=0-3;yo=  1)  +  He  collisions  in  the  cenler-of-mass  collision  energy  {Ej)  range  of  0.0-2. 0  eV. 
The  present  TWQS  approach  accounts  for  Coriolis  coupling  and  uses  the  ah  initio  potential  energy 
surface  of  Palmieri  et  al.  [Mol.  Phys.  98,  1839  (2000)].  For  a  fixed  total  angular  momentum  J,  the 
energy  dependence  of  reaction  probabilities  exhibits  quantum  resonance  structure.  The  resonances 
are  more  pronounced  for  low  J  values  and  for  the  HeH'^  +  D  channel  than  for  the  HeD'^+H  channel 
and  are  particularly  prominent  near  threshold.  The  quantum  effects  are  no  longer  discemable  in  the 
integral  cross  sections,  which  compare  closely  to  quasiclassical  trajectory  calculations  conducted  on 
the  same  potential  energy  surface.  The  integral  cross  sections  also  compare  well  to  recent 
state-selected  experimental  values  over  the  same  reactant  and  translational  energy  range.  Classical 
impulsive  dynamics  and  steric  arguments  can  account  for  the  significant  isotope  effect  in  favor  of 
the  deuteron  transfer  channel  observed  for  HD'^(t;<3)  and  low  translational  energies.  At  higher 
reactant  energies,  angular  momentum  constraints  favor  the  proton-transfer  channel,  and  isotopic 
differences  in  the  integral  cross  sections  are  no  longer  significant.  The  integral  cross  sections  as  well 
as  the  J  dependence  of  partial  cross  sections  exhibit  a  significant  alignment  effect  in  favor  of 
collisions  with  the  rotational  angular  momentum  vector  perpendicular  to  the  Jacobi  R 

coordinate.  This  effect  is  most  pronounced  for  the  proton-transfer  channel  at  low  vibrational  and 
translational  energies.  ©  2007  American  Institute  of  Physics.  [DOI:  10. 1063/1. 2800009] 


I.  INTRODUCTION 

The  endothermic  HD'^  +  He  reaction  to  form  the  isotopo- 
logues,  HeH^  and  HeD^ 


HeH*+D, 

A// =  0.8 1 7  eV, 

(1) 

HeD*  +  H, 

A// =  0.777  eV, 

(2) 

is  among  the  most  fundamental  chemical  reaction  systems.  It 
is  a  textbook  example  for  the  study  of  the  translational  and 
internal  energy  dependences  of  chemical  reactivity.  A  num¬ 
ber  of  experimental  and  theoretical  studies  have  investigated 
isotopic  branching  in  the  HD'^-l-He  reactive  collision 
system.*"*®  Significant  discrepancies  exist  between  theory 
and  state-selected  studies  by  Turner  et  al.,^  in  particular,  with 
regard  to  the  vibrational  and  translational  energy  depen¬ 
dences  of  integral  reaction  cross  section  branching  ratios, 
r(t;)  =  crj^(HeH^)/o‘j^(HeD^),  for  i;=0-4.  More  recently, 
Tang  et  al.^^  reported  absolute  integral  proton-transfer  and 
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deuteron  transfer  cross  sections  for  state-selected  reactants, 
HD'^(t;,y=  O  +  He,  using  the  vacuum  ultraviolet  (VUV) 
photoionization-guided  ion  beam  (VUV-PI-GIB)  technique 
and  the  pulsed  field  ionization-photoeicctron-sccondary  ion 
coincidence  scheme.  These  measurements,  which  signifi¬ 
cantly  extended  the  range  of  reactant  vibrational  and  transla¬ 
tional  energies  for  which  absolute  integral  cross  sections  are 
available  for  channels  (1)  and  (2),  provided  a  much  im¬ 
proved  comparison  of  determined  branching  ratios  at  a  trans¬ 
lational  energy  of  I  eV  with  both  quasiclassical  trajectory 
(QCT)  as  well  as  time-deprendent  wave  packet  quantum  scat¬ 
tering  (TWQS)  calculations  by  Tiwari  et  al.^  The  re.sults  by 
Tang  et  al.,  therefore,  suggest  that  the  measurements  by 
Turner  et  al.,  whose  branching  ratios  are  approximately  a 
factor  of  2  higher  than  tho.se  of  Tang  et  al.  and  Tiwari  et  al. 
for  t;  =  3  and  4,  are  inaccurate.  However,  the  new  experimen¬ 
tal  integral  cross  sections  were  in  much  better  agreement 
with  QCT  calculations  than  the  TWQS  results.  Additional 
quantum  studies  including  high  J  are,  therefore,  warranted  to 
resolve  the  discrepancies. 

The  Ht+He  collision  system  has  been  the  subject  of 
extensive  theoretical  work.^"®**^”^^  Of  particular  interest  has 
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been  the  observation  of  prominent  resonant  structure  in  the 
energy  dependence  of  reaction  probabilities  in  quantum  scat¬ 
tering  experiments.^’*^’^^  These  calculations  suggest  that 
the  QCT  method  may  not  be  appropriate  in  elucidating  the 
dynamics  of  this  fundamental  reaction  system.  A  recent  ac- 
curate  potential  energy  surface  calculated  by  Palmieri  ct  al: 
has  provided  a  new  impetus  to  theoretical  work  on  both  Ht 
-f-He  reactivity  as  well  as  the  isotopic  reactive  channels. 
Chu  et  used  the  Palmieri  surface  to  conduct  time- 
dependent  wave  packet  quantum  scattering  calculations  in¬ 
cluding  Coriolis  coupling  (TWQS-CC)  to  determine  absolute 
integral  Ht(t^=0-2,4,6,y=  l)-f  He  proton-transfer  cross  sec¬ 
tions  over  an  unprecedented  translational  energy  range,  Ej 
=  0-2.4  eV.  Their  results  demonstrated  that  the  inclusion  of 
the  Coriolis  coupling  is  e.ssential  for  obtaining  satisfactory 
agreement  between  theoretical  and  experimental  cross  sec¬ 
tions  for  HeH'^-nH  formation  from  the  Ht+He  reaction. 
Overall,  the  TWQS-CC  results  compared  well  to  the  most 
recent  state-selected  experiments  by  Tang  ct  al.^^  This  work 
also  demon.strated  that  despite  the  prominence  of  resonances 
in  y-selected  reaction  probabilities,  the  TWQS-CC  integral 
cross  sections  compared  very  well  to  QCT  calculations  using 
the  same  potential  energy  surface. 

In  this  paper,  wc  present  a  TWQS-CC  study  of  the  iso¬ 
topic  reaction  channels  in  HD^(i;  =  0-3  ;y=  l)-»-He  collisions 
for  the  collision  energy  range  of  0.0-2. 0  eV.  A  detailed  com¬ 
parison  is  provided  between  the  pre.senl  TWQS-CC  results, 
recent  experimental  work  reported  by  Tang  et  r//.,**  and  QCT 
calculations  for  similar  reactant  energies. 


II.  THEORETICAL  METHOD 

The  prc.scnt  TWQS-CC  study  of  reactions  (1)  and  (2)  is 
similar  to  the  recent  theoretical  investigation  of  the  proton- 
transfer  reaction  Ht(  t^,y)-t-He— HeH^-i-H  reported  by  Chu 
et  The  time-dependent  wave  packet  method  has  been 
widely  u.sed  in  dynamical  studies  for  the  elucidation  of  the 
mechanism  of  simple  reactive  systems."^*  This  methexi  nu¬ 
merically  solves  the  time-dependent  Schrodinger  equation  of 
the  reactive  system  in  conjunction  with  the  split-operator 
propagator  scheme. 


In  the  reactant  Jacobi  coordinates,  the  Hamiltonian  of 
the  HeHD"^  system  can  he  expressed  as 


f,2  ^  f 


V(R.r)-\-h(r). 


(3) 


where  R  is  the  distance  between  the  center  of  mass  of  HD^ 
and  He,  r  is  the  HD'*‘  bond  length,  is  the  reduced  mass  of 
He  with  respect  to  HD"^,  is  the  HD'*‘  reduced  mass,*^  J  and 
j  are  the  respective  total  angular  momentum  of  the  HD"^ 
-»-He  reaction  system  and  the  rotational  angular  momentum 
of  HD"^,  V(R,r)  is  the  |)otential  energy  surface,  as  deter¬ 
mined  by  Palmieri  et  and  h{r)  is  the  diatomic  reference 
Hamiltonian, 


/Kr)  =  - 


tr  iP- 
2m,  <)r 


(4) 


where  V{r)  is  the  diatomic  reference  potential. 

The  time-dependent  wave  function  is  expanded  in  terms 
of  the  translational  basis  the  vibrational  basis 

and  the  hody-lixed  total  angular  momentum  basis 
where  M  and  K  are  the  projection  of  J  on  the 
space- fixed  and  body- fixed  r.  axis,  respectively,  and  e  pro¬ 
vides  the  total  parity.  In  the  coupled-state  or  centrifugal  sud¬ 
den  approximation,*'*^*'^  the  matrix  obtained  by  applying  the 
centrifugal  potential  operator,  the  second  term  of  the  Hamil¬ 
tonian  [Eq.  (3)],  to  the  rotational  basis  functions  Y\^J^{R,r)  is 

'  ^  44  4S 

assumed  diagonal,  whereas  in  the  pre.sent  CC  calculation, 
different  K  states  are  allowed  to  couple,  leading  to  off- 
diagonal  elements. 

The  initial  wave  function  is  propagated  by  a  split 
operator  scheme,**^  and  the  ./,  y'o,  A(),  and  v  selected  reaction 
probability  the  yj),  k^),  and  u- selected  integral 

reaction  cross  section  rovihra- 

tionally  state-selected  integral  reaction  cross  sections 
calculated  from*** 


004 


0  0  0  5  TO  1  5  0.0  0  5  1  0  1  5  0  0 

E,  (eV) 


0  03 

002 

001 

0,00 
0  12 

000 


(b)  HeH+.  ko=-1  (c)  ko=0 


J  (a)HeH^, 

JO 


(d) 


(e) 


(f) 


FIG.  I.  (Color  onlinci  TWQS-CC  re- 
aclion  channel  probabililics. 
(u=0,7()=l)  for  y=  1  (black  curves), 
10  (red  curves).  20  (green  curves),  30 
(pink  curves),  40  (blue  curves),  ami  50 
(dark  green  curves)  in  the  range  o\ 
().()-2.()eV  for  Ao=().±l. 
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FIG.  2.  (Color)  TWQS-CC  rcaclion 
channel  probabiliiies. 

{v~  \  Jq^\)  for  y=l  (black  curves), 
10  (red  cur\es),  20  (green  curves).  30 
(pink  curves).  40  (blue  curves),  and  50 
(dark  green  curves)  in  ihe  range  of 
0.0-2.0  eV  for  A„=0.±1. 


^  Im[<^4(£r)|<5(r-  ro)|^i(£r))]-  (5) 

(2i  +  .  (6) 

^  J 

=  ^777^  .  (7) 

JO  k^^ 

where  (piEj)  corresponding  time-independent  part  of 

the  final  wave  function,  ^o  is  the  projection  of  the  initial 
rotational  angular  momentum  Jq  on  the  body-fixed  z  axis, 
and  k={2EjiJir^fry^^  is  the  wave  number  corresponding  to 
the  initial  state  at  a  fixed  collision  energy,  Ej. 

The  present  theoretical  work  is  also  complimented  with 
QCT  calculations  using  the  same  potential  energy  surface. 
The  present  QCT  approach  has  been  described 
previou.sly.*^*"*^  While  most  of  the  QCT  results  have  been 
reported  already  by  Tang  et  ai^  additional  calculations  were 
conducted  for  those  comparisons  where  higher  statistics  were 
necessary.  Up  to  90  000  trajectories  were  computed  in  those 
ca.ses. 


III.  RESULTS 
A.  Reaction  probabilities 

Figures  1(a)- 1(c)  [1(d)- 1(f)]  depict  the  translational  en¬ 
ergy  dependencies  between  0  and  2.0  eV  of  the  TWQS-CC 
reaction  probabilities  for  the  formation  of  HeH"^ 

[HeD^J  calculated  with  HD"^  in  the  respective  (t;=0; 

^0=0.  ±1)  states  at  J=  I  (black  curves).  10  (red  curves),  20 
(green  curves),  30  (pink  curves),  40  (blue  curves),  and  50 
(dark  green  curve.s).  The  corresponding  translational  energy 
dependencies  of  TWQS  CC  reaction  probabilities  for  t;=l, 
2,  and  3  calculated  for  the  same  J  values  are  compared  in 
Figs.  2-4,  respectively. 

Sharp  resonance  structure  is  observed,  particularly  for 
J=  I  collisions  near  threshold.  As  J  is  increased,  the  reso¬ 
nance  structure  becomes  less  evident.  Reactive  resonances 
are  more  pronounced  for  HeH"^  than  for  HeD"^.  The  general 
trends  of  the  probabilities  with  energy  are  also  quite  different 
for  the  two  isotopic  channels.  The  reaction  probabilities  for 
HeD"^  increase  gradually  as  a  function  of  Ej  except  near  the 
probability  maximum  of  7^20,  whereas  the  probabilities  for 
HeH"^  generally  have  a  sharper  onset  at  threshold,  above 
which  the  resonances  are  the  most  pronounced.  For  a  fixed 
reactant  vibrational  level  and  a  fixed  J  value,  the  behaviors 
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FIG  3.  (Color)  TWQS-CC  rcaclion 
channel  probabiliiies. 

(i;=2,)o=l)  for  7=1  (black  curves). 
10  (red  curves).  20  (green  curves),  .30 
(pink  curves).  40  (blue  curvc.s).  and  .50 
(dark  green  curve.s)  in  ihe  range  of 
0.0-2.0  eV  for  ito=0.±F 
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MG.  4.  (Color)  TWQS  CC  rciiclion 
channel  probabiluies. 

(u  =  3.7o=I)  (black  curves). 

10  (red  curves),  20  (green  curves).  30 
(pink  curves),  40  (blue  curves),  and  50 
(dark  green  curves)  in  ihe  Ej  range  ot 
00-2,0  cV  for  k^,={).±\. 


of  the  reaction  probabilities  for  k^)=  1  and  A()=-l  are  similar, 
but  are  substantially  different  from  the  behavior  for  kQ=0 
reactants.  The  probability  curves  of  both  channels  calculated 
for  ko=±\  and  i;=()-3  are  shown  to  decrease  as  J  is  in¬ 
creased,  whereas  the  Ao=0  probability  curves  of  HeD'*’  for 
i;=0-3  and  HeH"^  for  i;=l-3  have  considerable  contribu¬ 
tions  from  high  J  collisions.  We  note  that  since  the  reaction 
cross  section  is  assexiated  with  a  reaction  probability  scaled 
by  (27+1),  the  contribution  to  the  reaction  cross  sections 
from  higher  7  values  is  larger  than  appears  from  the  reaction 
probabilities.  Furthermore,  for  all  k^^  and  7  values,  the  Ej 


increased.  This  trend  is  more  pronounced  as  the  reactant  vi¬ 
brational  energy  increases. 

A  more  comprehensive  three  dimensional  (3D)  represen¬ 
tation  of  .Ef)  is  shown  for  i;=0  and  1  reactant  ions  in 

Figs.  5  and  6,  respectively.  The  k[)=-\  states  are  not  shown 
because  they  are  very  similar  to  the  corresponding  plots  for 
A'()=l  states.  The  3D  graphs  rcalfirm  the  above  observation 


that  reactive  resonances  are  significantly  sharper  and  more 
pronounced  for  the  HeH'*’  channel.  The  resonance  structure 
diminishes  as  i;  is  increased  from  0  to  3.  Corresponding  plots 
for  u  =  2  and  i;  =  3  can  be  obtained  from  the  authors. 

B,  J  dependence  of  partial  cross  sections 

The  proton  or  deuteron  transfer  TWQS-CC  partial  cross 
sections  for  a  given  reactant  state  HD'^(v 
f/y^^()t^^(7,£V)  =  (7r/A“)(27+  l)/^^^^^^,(£Y),  depend  on  both  7  and 
Ef.  Figure  7  exhibits  the  7  dependencies  of  the  respective 
partial  cross  sections  for  i;  =  ()-3;  yo=l,  A()=0,  ±1  at  Ep 
=  I  eV  for  the  two  isotopic  channels.  These  partial  cross  sec¬ 
tion  plots  show  that  the  theoretical  partial  cross  sections  for 
k^^=  1  and  —1  arc  cssciuially  identical  and  that  the  'fWQS-CC' 
cross  sections  at  Ej=  I  eV  vanish  above  7^70.  A  signiheant 
preference  for  Ao=()  reactants  is  observed  in  the  HeH^  +  D 
channel  for  all  investigated  reactant  states  and  for  i;=()  reac¬ 
tants  in  the  HeD‘*’  +  H  channel.  For  i;>0  reactants,  the 
HeD'*’+H  partial  cross  sections  at  1  eV  are  nearly  indepen- 


FIG.  5.  (Color)  The  reaclion  probabililies  (u  =  0.7()=  I  ./(o^O.  F-IG.  6.  (Color)  Tbe  reaclion  probabililies  iv=  \  .j=  \  .k=i). 

+  I )  for  HeH^-i-D  and  HeD*+  H  formation  plotted  a.s  a  function  of  J  and  ty  +  1 )  for  HeH*  +  D  and  HeD*  +  H  formation  plotted  a,s  a  function  of  J  and  Er 
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FIG.  7.  J  dependencies  of  ihe 
TWQS-CC  partial  cross  sections, 
(jo=l./-'r=I  eV)  for  t; 
=  0-3,  and  Ao  =  0,  ±  1  for  the  fomiation 
of  HeH*  (solid  curves)  and  HeD* 
(dashed  curves). 


dent  of  As  expected,  given  the  endothermicity  of  the  re¬ 
action,  the  TWQS-CC  partial  cross  sections  for  HeH"^  and 
HeD^  increa.se  .significantly  as  v  increases.  There  is  a  gradual 
shift  of  the  J  distribution  to  higher  values  as  v  increases. 
Furthermore,  except  for  i;=0  reactants,  there  is  a  clear  pref¬ 
erence  for  HeD^  formation  at  low  J  values,  whereas  the  for¬ 
mation  of  HeH"^  is  favored  at  high  J  values.  The  cross-over  J 
value  above  which  the  HeH'*‘  channel  is  preferred  occurs  be¬ 
tween  35  and  40  and  does  not  show  a  clear  dependence  on 
i;  for  v=  1  -3. 

Figure  8  compares  the  /:o-averaged  TWQS-CC  partial 
cross  sections  with  QCT  calculations  based  on  the  same  po- 
tential  energy  surface  by  Palmieri  et  al.  The  previous  com¬ 
parison  between  experimental  and  QCT  integral  cross 
.sections**  has  demonstrated  that  the  agreement  is  most  favor¬ 
able  if  the  zero-point  energy  in  reaction  products  is  disre¬ 
garded.  In  this  work,  we  consequently  only  compare  with 
QCT  calculations  that  do  not  account  for  product  zero-point 
energy.  The  absolute  (JCT  values,  shown  as  solid  and  dashed 
lines  for  HeH"^  and  HeD"^,  respectively,  are  derived  from 
fiL,  where  L=fji^gb  is  the  classical  correspondence 
to  1)A^  assuming  y=|l+j|=*®/,  and  g  and  h  are  the 

relative  velocities  at  infinite  distance  and  impact  parameter, 
respectively.  There  is  generally  a  good  agreement  between 
the  classical  and  TWQS-CC  partial  cross  sections,  both  in 
magnitude  and  J  dependence.  The  classical  distributions  ex¬ 
hibit  a  minor  shift  toward  higher  angular  momentum  colli¬ 
sions  in  both  channels.  Probably  the  most  significant  differ¬ 
ences  between  TWQS-CC  and  QCT  are  seen  in  the  HeH'*’ 
channel  for  t;=0  reactants,  where  unlike  QCT,  the 
TWQS-CC  results  exhibit  bi modality  with  a  tail  toward 
high  y. 


C.  Rovibrationally  state-selected  integral  cross 
sections 

The  TWQS-CC  integral  cross  .sections, 

=  0-3\jQ=\  \k^^=0,±\),  for  both  isotopic  channels  calcu¬ 
lated  for  £7^= 0.0-2. 0  eV  are  plotted  in  I^g.  9.  The  solid 
curves  shown  in  these  figures  are  the  ^Q-averaged,  rovibra¬ 
tionally  St  ate- selected  integral  cross  sections,  (Tj^^^X^ry 
seen  that  the  A'o=0  integral  cross  sections  for  HeD”^  and 
HeH"^  are  larger  than  the  corresponding  ^0  =  ±  1  cro.ss  sec¬ 
tions  over  the  whole  Ej  range  except  the  HeD‘*‘  values  for 
v  =  2  and  3,  where  the  cross  sections  are  observed  to  be  es¬ 
sentially  independent  of  Aq.  All  the  cro.ss  section  curves  ex¬ 
hibit  relatively  sharp  Ej  thresholds.  The  pronounced  reso¬ 
nance  structure  observed  in  the  probability  curves  is  greatly 
suppressed  in  the  TWQS-CC  integral  cross  sections  due  to  J 
and  A'o  averaging.  With  the  exception  of  the  t;=0  integral 
cross  sections  of  the  HeD"^  channel  [Fig.  9(a)],  which  still 
increase  at  Ej=2  eV,  the  HeD'^  integral  cro.ss  .sections  for 
reactant  states  t;=l-3  [Figs.  9(b)-9(d)]  and  those  of  HeH"^ 
for  t;=0-3  [Figs.  9(aa)-9(dd)]  are  found  to  decline  at  ener¬ 
gies  above  I  eV.  The  decline  is  found  to  be  more  rapid  in  the 
case  of  HeH'*'  formation.  Also,  the  integral  reaction  cross 
sections  for  the  HeH"^  isotopic  channel  display  sharper 
maxima.  Interestingly,  the  HeD'^  cross  sections  for  t;=l-3 
reactants  [Figs.  9(b)-9(d)]  manifest  some  structure  near 
threshold. 

Figure  9  also  compares  the  Ag-averaged  TWQS-CC  cross 
sections  with  QCT  calculations  using  the  surface  by  Palmieri 
et  Except  close  to  threshold,  the  QCT  predictions  are 
generally  in  very  good  agreement  with  the  TWQS-CC  re¬ 
sults.  The  discrepancy  at  threshold  is  expected  because  zero- 
point  energy  of  product  molecular  ions  is  not  taken  into 
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MG,  8.  /to-avcniged  TWQS-CC  partial  cross  scclion,s  lor  Ej=  I  eV. 

70=1.  and  (;  =  0-3  (filled  circlc,s:  HelT  +  D:  open  circles:  HeD*  +  H).  The 
solid  and  broken  lines  are  the  corresponding  C)CT  re.sults, 

account  in  the  QCT  calculations.  Consequently,  the  QCT 
thresholds  are  not  just  lower  than  the  exact  thresholds,  but 
are  identical  for  the  two  isotopic  channels. 


D.  Comparison  with  experimental  cross  sections 

A  detailed  comparison  of  existing  experimental  and  the¬ 
oretical  cross  .sections  for  reactions  ( 1 )  and  (2)  was  recently 
reported  by  Tang  ei  r//.  **  In  Fig.  10  we  compare  the  VUV- 
PI-GIB  measurements  of  integral  reaction  cross  sections  re¬ 
ported  by  Tang  et  ai  to  the  present  TWQS-CC  re.sults.  The 
theoretical  cross  section  energy  dependence  has  been  convo¬ 
luted  by  the  experimental  energy  broadening  functions,  in¬ 
cluding  both  thermal  target  gas  broadening^^  and  the  ion  en¬ 
ergy  resolution  of  the  experiment  (~0.3  eV).  Any  structure 
observed  in  the  TWQS-CC  cross  section  energy  dependen¬ 
cies,  particularly  the  shoulder  observed  in  the  HeD^  onsets, 
is  no  longer  apparent  following  the  convolution.  Above 
1  eV,  the  agreement  between  experimental  and  theoretical 
cross  sections  is  excellent.  At  lower  energies,  the  experimen¬ 
tal  values  are  slightly  higher  than  the  theoretical  predictions, 
particularly  near  the  cross  section  maxima  for  t;=l  and  2. 
For  all  reactant  states,  the  energy  regions  of  strong  cross 
section  growth  with  energy  appear  at  slightly  higher  energies 


0.0  0.5  1  0  1  5  0  0  0  5  1.0  1.5  2,0 


E,  (eV) 

MG.  9.  The  TWQS-CC  prolon  and  deulcron  Iransfcr  integral  cross  seclions 
XT with  7,)=  1  and  (dashed  line),  ±1  (dotted  lines),  and 

/t'o-ave raged  (solid  line)  in  the  Ej  range  of  0.0  2.0  eV  for  t;  =  0  3.  The 
TWQS-CC  results  are  compared  with  QCT  ero.ss  scelions. 

in  the  case  of  the  theoretical  results.  This  shift  in  apparent 
cross  section  onset  is  more  pronounced  in  the  case  of  the 
deuteron  transfer  channel. 

The  source  of  the  discrepancy  near  threshold  could  be 
caused  by  translational  energy  calibration  error  in  the  experi 
ments,  or  by  differences  in  the  steepness  of  the  experimental 
and  theoretical  threshold  onsets,  or  by  inaccuracies  in  the 
potential  energy  surface.  Tang  ei  alV  found  that  the  ob 
served  cross  .section  energy  dependences  for  the  four  lowest 
reactant  vibrational  .states  were  consistent  with  modified 
line-of-center***^***^  threshold  functions  based  on  the  known 
thermochemical  thresholds  of  the  two  isotopic  channels.  The 
deconvoluted  cross  .sections,  however,  were  found  to  have 
highly  vertical  onsets  at  threshold.  A  difference  in  threshold 
function  steepness  could  also  be  related  to  the  fact  that  the 
VUV-PI-GIB  measurements  by  Tang  et  ai  were  not  rotation- 
ally  state  selected,  unlike  the  TWQS-CC  calculations  which 
were  conducted  for  7o=l.  This  is  further  discussed  in  the 
next  section. 

Table  I  compares  the  branching  ratios  FCu)  derived  from 
the  present  calculations  at  I  eV  to  those  obtained  from  the 
present  QCT  calculations  and  other  theoretical  work,  as  well 
as  experimental  values.  The  TWQS-CC  results  are  in  gotxl 
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FIG.  10.  Compari.son  between  the  energy  dependence  of  experimental  (filled  and  open  circles)  and  TWQS-CC  (.solid  and  da.shed  linc.s)  integral  reaction  cross 
sections.  The  TWQS-CC  curves  have  been  convoluted  by  the  experimental  energy  broadening  function,  including  primary  ion  energy  distribution  and  target 
gas  thermal  motion. 


agreement  with  the  experimental  values  apart  from  slightly 
higher  values  at  i;=0  and  2.  Except  at  t;=0,  the  present  QCT 
results  have  higher  ratios  than  those  produced  by  TWQS-CC. 
The  TWQS  branching  ratios  by  Tiwari  et  arc  signifi¬ 
cantly  different  from  the  experimental  ones  for  v<2. 

IV.  DISCUSSION 

The  present  work  confirms  the  accuracy  of  the 
TWQS-CC  method  in  predicting  accurate  integral  cross  sec¬ 
tions  for  the  textbook  H2(HD^)  +  He  reaction  systems  at  total 
energies  as  high  as  2.8  eV  using  the  potential  energy  surface 
of  Palmieri  et  al.^^  The  isotopic  branching  in  the  present 
study  is  a  particularly  sensitive  probe  of  the  dynamics.  As 


seen  in  Table  I  and  Fig.  10,  the  calculations  are  successful  in 
reproducing  the  vibrationally  state-selected  experimental  in¬ 
tegral  cross  sections  and  respective  isotopic  ratios  by  Tang  et 
al.^^  at  translational  energies  above  1  eV.  Closer  to  thresh¬ 
old,  the  calculations  appear  to  underpredict  the  i;  =  2  and  3 
HeD^  channel  cross  sections;  however,  the  theoretical  results 
are  still  within  the  stated  experimental  accuracy.  The  latter  is 
reduced  near  threshold  due  to  the  relatively  broad  experi¬ 
mental  translational  energy  distribution  in  a  region  of  strong 
cross  section  energy  dependence.  As  mentioned  in  the  previ¬ 
ous  section,  Tang  et  t//. "  deconvoluted  the  experimental 
cross  sections  to  derive  modified  line-of-center  threshold 
functions  that  have  significantly  steeper  onsets  than  predicted 
by  the  present  theory.  Assuming  that  no  translational  calibra- 


TABLE  I.  Comparison  of  the  experimental  and  theoretical  values  for  the  branching  ratios  \'{v) 
=<rj,(HeH^)/<ry(HeD^)(i;=0-3)  obtained  at  £‘7-=1.0eV,  where  and  <ry(HeD^)  are  integral  cro.ss 

section.s  for  the  formation  of  HeH^  and  HeD^  respectively. 


state 

Theory 

Experiment* 

TWQS-CC^ 

TWQS^ 

QCT* 

ocr* 

TSlf 

o 

II 

0.39 

0.98 

0.54 

1.06 

0.35±0.04 

0.34 

0.75 

i;=I 

0.51 

0.72 

0.47 

0.23 

0.69±0.06 

0.92 

0.33 

i;  =  2 

0.57 

1.04 

0.77 

0.52 

0.91  ±0.04 

0.92 

0.71 

II 

0.70 

1.38 

0.83 

0.83 

1.11±0.04 

1.20 

1.20 

“Left:  Ref.  1 1;  right:  Ref.  4. 
^Present. 


‘^Reference  9. 
‘^Reference  5. 
^Reference  6. 
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FIG.  11.  Schematic  illusiralion  of  ^=1  projections.  A.'q=0,±1,  where  R  is 
Ihc  quantization  axis. 


tion  error  is  the  cause  of  the  discrepancy,  it  is  possible  that 
the  population  of  reactant  rotational  slates  ^  1  in  the  ex¬ 
periments  and  a  particular  sensitivity  to  reactant  rotational 
energy  is  responsible  for  the  steeper  experimental  onset.  In- 
deed,  Kalyanaraman  et  cil  discovered  a  strong  y’o  depen¬ 
dence  of  i=0  reaction  probabilities  in  quantum  scattering 
calculations.  They  found  that  the  HeH'^-i-D  channel  70  depen¬ 
dence  was  irregular  with  t;,  while  the  HeD'^  +  H  channel  ex¬ 
hibited  a  preference  for  even  y’o  quantum  numbers  for  all  v. 
with  reactants  exhibiting  the  lowest  probabilities.  Sub¬ 
sequently,  Tiwari  and  Saihyamurthy*^  calculated  the  rota¬ 
tional  energy  dependence  of  integral  cross  sections  for  1;=  I 
reactants  using  a  three  dimensional  TWQS  approach.  Their 
results  show  that  in  the  case  of  the  proton-transfer  channel, 
7o=  1  cross  sections  are  slightly  lower  than  those  for  jQ=0, 
but  significantly  higher  than  7o=2  and  3,  while  the  deuleron 
transfer  cross  sections  of  1  reactants  are  substantially 
lower  than  those  of  70=0,  2,  and  3.  This  observation  would 
be  consistent  with  observation  that  the  theoretical  70=  1  deu- 
teron  transfer  channel  cross  sections  exhibit  larger  differ¬ 
ences  with  the  experimental  cross  sections  than  the  proton- 
transfer  results  near  threshold. 

Additional  evidence  for  reactant  rotational  sensitivity  is 
provided  by  the  significant  Aq  dependence  of  the  HeH'^-i-D 
channel  integral  state-selected  cross  sections  at  all  investi¬ 
gated  reactant  vibrational  levels  and  for  i;=0  and  1  reactants 
for  the  HeD"^  channel  (Fig.  9),  Similar  observations  were 
made  by  Tiwari  and  Salhyamurty  for  i;=l  reactants.*^  The 
favoring  of  Ao=0  reactants  signifies  a  preferred  reactant 
alignment  with  the  rotational  axis  perpendicular  to  the 
Jacobi  R  coordinate,  represented  by  the  line  between  the  He 
target  and  the  center  of  mass.  This  is  shown  schemati¬ 
cally  in  Fig.  1 1  for  collisions  at  the  point  of  approach,  R 
=»0.7  A,  which  is  also  the  equilibrium  distance  between  the 
H  atom  and  the  center  of  mass  of  HD"^.  The  preference  of 
Ao=0  is  consistent  with  the  minimum  energy  reaction  path 
proceeding  through  a  col  linear  geometry. 

In  Fig.  7,  it  is  seen  that  at  I  eV,  the  i;=0,  Ao=0  partial 
cross  sections  of  both  isotopic  channels  are  largest  near  J 
=  20  which  corresponds  to  an  impact  parameter  of  -^0.7  A. 
For  Ao=±l,  the  distributions  peak  near  7=15.  At  all  Aq,  v 
=0  reactants  preferably  form  HeD'^-i-H  products.  This  is  con¬ 
sistent  with  the  necessity  to  transfer  -^80%  of  the  transla¬ 


tional  energy  into  internal  energy  through  small  impact  pa¬ 
rameter  collisions,  whereby  the  impact  parameters  have  to  be 
slightly  smaller  for  the  Ao=±l  rotational  orientations  for 
H/D  atoms  to  undergo  sufficiently  close  interaction  with  the 
He  target  atom. 

The  excellent  agreement  between  QCT  and  the  present 
quantum  calculations  suggests  that  the  minor  difference  in 
isotopic  thresholds  only  plays  a  role  in  the  preference  for 
HeD'^-i-He  products  at  energies  above  threshold  comparable 
to  the  isotopic  difference.  The  favoring  of  the  deuteron  trans 
fer  channel  can  be  much  more  readily  explained  using  en¬ 
ergy  transfer  arguments  associated  with  impulsive  or 
spectator- St  ripping- type  dynamics.  Tang  et  interpreted 
recoil  velocity  distributions  observed  for  i;  =  3  reactants  at 
£^=0.6  eV,  as  evidence  for  an  impulsive  or  spectator- 
stripping  mechanism.  In  a  purely  impulsive  mechanism  in 
which  only  the  transferred  proton/de uteron  collides  with  the 
He  target  atom,  47%  of  the  translational  energy  is  transferred 
into  the  product  ion  in  the  case  of  the  proton-transfer  chan¬ 
nel,  while  78%  is  transferred  in  the  case  of  the  deuteron 
transfer  channel.  Thus,  for  endothermic  prex'esses,  a  direct 
mechanism  would  strongly  favor  the  HeD'^+H  products.  In 
the  case  of  the  high-threshold  i;=0  process,  additional  favor¬ 
ing  of  the  deuteron  transfer  channel  may  be  attributable  to 
the  on  average  closer  He-D  distance  at  the  turning  p<nnt  of 
the  collision  given  the  tighter  gyrating  trajectory  of  the  D 
atom  due  to  its  closer  vicinity  to  the  center  of  mass  of  the 
rotating  HD"*^  molecule. 

For  i;=l  reactants,  the  HeH"^  channel  is  preferred  by 
higher  J  collisions  with  still  a  strong  preference  for  the  A^ 
=0  rotational  alignment,  while  the  alignment  effect  is  less 
pronounced  for  the  HeD"^  channel  (see  Fig.  7).  For  the  latter, 
the  partial  cross  sections  still  peak  around  7  =  20,  with  the 
A()=0  distribution  at  slightly  higher  7  than  the  Ao=±  1  distri¬ 
butions.  As  pointed  out  by  Tang  et  «/.,"  the  significant  dif¬ 
ference  in  the  opacity  functions  of  the  two  isotopic  channels 
can  be  rationalized  by  the  significant  difference  in  the  re¬ 
duced  mass  of  the  two  product  channels  coupled  with  the 
steric  argument  of  the  longer  reach  of  the  H  atom  with  re¬ 
spect  to  the  HD"^  center  of  mass.  While  the  HeH'^  +  D  channel 
has  a  reduced  mass  of  1.43  amu  that  is  similar  to  that  of  the 
reactants  (1.71  amu),  the  reduced  mass  of  HeD'^  +  H  products 
is  significantly  lower  (0.86  amu)  than  those  of  the  reaetanis, 
and  thus,  substantial  angular  momentum  conservation  con¬ 
straints  arise  at  large  7  for  the  deuleron  transfer  channel. 
Interestingly,  Light  and  Lin  have  predicted  this  large  isotope 
effect  in  the  opacity  functions  for  this  reaction  due  to  angular 
momentum  conservation  arguments  more  than  40  years  ago.^ 

The  large  preference  for  Ao=0  reactants  in  the  HeH'^ 
channel  can  be  explained  by  the  higher  probability  of  trans 
ferring  sufficient  energy  in  the  large  impact  parameter  colli¬ 
sions  that  appear  to  be  dynamically  preferred  by  this  chan¬ 
nel.  This  is  also  seen  to  be  the  case  for  u  =  2  and  3  reactants. 
However,  in  the  case  of  HcD'^  +  H  products,  there  is  no 
longer  a  significant  difference  between  the  Ao-selected  inie 
gral  cross  sections,  and  only  a  minor  shift  in  the  7  disiribu 
lions  is  observed.  This  is  consistent  with  a  7  distribution  that 
is  primarily  dictated  by  angular  momentum  constraints  ver¬ 
sus  energy  transfer  requirements.  The  latter  are  not  as  sirin- 
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gent  for  the  deutcron  transfer  channel  due  to  the  higher  effi¬ 
ciency  in  impulsive  dynamics.  The  present  classical 
rationalization  is  further  justified  by  the  close  resemblance  of 
the  k^-averiiged  partial  cross  section  distributions  and  corre¬ 
sponding  QCT  results  shown  in  Fig.  8.  Qualitatively,  the 
partial  cross  sections  calculated  by  Tiwari  et  ai  '  using 
TWQS  with  both  the  Palmieri  et  and 

McLaughlin-Thompson-Joseph-Sathyamurthy  ’  potential 
energy  surfaces  are  similar  to  the  present  distributions;  how¬ 
ever,  the  present  calculations  including  Coriolis  coupling  re¬ 
semble  the  QCT  predictions  more  closely  at  £7-=  1  eV. 

Consistency  with  the  above  interpretations  is  also  found 
with  the  overall  behavior  of  the  translational  energy  depen¬ 
dence  of  the  /-dependent  reaction  probabilities  shown  in 
Figs.  1-6.  Three  types  of  dynamics  can  be  identified:  Sharp 
onsets  at  threshold;  dynamics  associated  with  a  steady  in¬ 
crease  followed  by  a  decrease  (if  extended  to  a  sufficiently 
high  in  cross  section  with  energy;  and  dynamics  associ¬ 
ated  with  resonance-type  structure.  Where  a  steady  increase 
in  cross  section  with  energy  is  ob.served,  this  behavior  turns 
on  at  lower  energies  in  the  case  of  the  HeD'*'-i-H  channel, 
which  is  consistent  with  impulsive  dynamics  and  the  higher 
efficiency  in  energy  transfer  for  this  channel.  Sharp  reso¬ 
nance  structure  due  to  quantum  effects  is  prominent  near 
thre.shold.  This  structure  is  particularly  evident  in  the  case  of 
the  HeH"^  channel  where  the  total  probabilities  are  signifi¬ 
cantly  lower.  Interestingly,  the  resonances  are  sharper  and 
more  intense  for  low  7,  kQ=±\  collisions,  suggesting  that 
they  are  less  likely  to  be  associated  with  the  attractive  well 
found  in  the  collinear  geometry.  The  resonances  are  more 
likely  attributable  to  trajectories  with  a  T-type  entrance  ge¬ 
ometry  that  are  associated  with  a  barrier.  In  this  geometry, 
the  interaction  of  He  is  more  likely  to  involve  both  H  and  D, 
contrary  to  more  impulsive  dynamics  involving  interaction 
between  the  He  atom  and  a  single  atom  of  the  HD"^  reactants 
in  the  collinear  geometry.  The  broader  resonances  ob¬ 
served  for  both  channels  at  higher  J  and  v  at  all  orientations 
are  .symptomatic  of  shorter-lived  intermediates.  As  already 
pointed  out  by  Light  and  Lin,^  statistical  density  of  state 
arguments  favor  the  deuteron  transfer  channel.  Thus,  the  de¬ 
cay  of  transition  state  resonances  would  occur  more  rapidly 
toward  deuteron  transfer  products,  explaining  the  broader 
resonances,  while  longer-lived  resonances  in  a  highly  indi¬ 
rect  mechanism  provide  more  efficient  energy  transfer  in  the 
proton-tran-sfer  channel.  For  the  latter,  differences  are  ob¬ 
served  between  kQ=-¥  1  and  -1  scattering  in  the  detailed  po¬ 
sitions  of  the  quantum  scattering  resonance  peaks,  suggest¬ 
ing  orientation  effects  at  this  level  of  detail. 

As  shown  in  Fig.  9,  resonances  do  not  appear  to  have  an 
important  effect  on  the  overall  reactivity  given  the  very  good 
agreement  between  the  quantum  and  QCT  inte¬ 

gral  cross  sections.  This  is  not  too  surprising  given  the  fact 
that  the  resonances  are  most  prominent  for  low  J  collisions, 
which  do  not  contribute  strongly  to  the  integral  cross  section. 
The  fact  that  the  QCT  calculations  that  do  not  take  zero-point 
energy  into  account  compare  so  well  to  the  present 
TWQS-CC  calculations  further  suggests  that  the  observed 
i.sotope  effects  shortly  above  threshold  have  a  purely  dy¬ 


namic  origin  and  that  the  reaction  threshold  difference, 
which  is  zero  for  the  QCT  calculations,  is  only  of  importance 
at  and  immediately  above  threshold. 

V.  CONCLUSIONS 

We  have  carried  out  a  3D  TWQS-CC  calculation  for 
the  rovibrationally  state-selected  isotopic  reactions  HD^ 
(i;=0-3;yo=l)  +  He->HeH^(HeD'')  +  D  (H)  in  the  Ej  range 
of  0.0-2. 0  eV  using  the  ab  initio  potential  energy  surface  of 
Palmieri  et  ai  Particular  attention  was  paid  to  the  rotational 
orientation  with  respect  to  the  body-fixed  c  axis.  The  results 
are  compared  with  QCT  calculations  on  the  same  potential 
energy  surface  and  the  state-selected  experiments  by  Tang  et 
ai^^  TWQS-CC  A-'o-averaged  integral  cross  sections  compare 
closely  to  QCT  predictions  over  most  of  the  investigated 
translational  energy  range,  with  only  important  discrepancies 
near  threshold.  The  Ao-averaged  partial  cross  sections 
at  1  eV  are  also  in  good  qualitative  agreement  with  QCT 
predictions,  suggesting  that  most  of  the  salient  dynamics  can 
be  explained  with  classical  arguments.  The  TWQS-CC  inte¬ 
gral  cross  sections  are  compared  to  state-selected  experi¬ 
ments  by  Tang  et  ai^^  following  convolution  of  the  quantum 
results  by  the  respective  experimental  energy  distribution 
function.  Above  threshold,  there  is  excellent  agreement  be¬ 
tween  the  experimental  and  theoretical  results.  The  experi¬ 
mental  energy  region  of  strong  cross  section  growth  with 
energy,  however,  is  slightly  shifted  to  lower  energies  than  the 
corresponding  convoluted  theoretical  values.  This  could  be 
attributed  to  a  lower  steepness  of  the  theoretical  cross  section 
onset,  possibly  attributable  to  important  rotational  effects 
given  the  fact  that  the  experiment  included  contributions 
from  thermally  populated  7*0  ^  *  reactant  states. 

TWQS-CC,  QCT,  and  experiment  demonstrate  a  signifi¬ 
cant  isotope  effect  for  i;<3  favoring  the  deuteron  transfer 
channel.  The  i.sotope  effect  is  primarily  the  result  of  compe¬ 
tition  between  impulsive  dynamics  favoring  the  deuteron 
transfer  channel  in  case  of  high-threshold  processes,  and  an¬ 
gular  momentum  constraints  favoring  proton-transfer  reac¬ 
tions  at  high  7.  Consequently,  indirect  dynamics,  including 
quantum  scattering  resonances,  play  an  important  role  in  the 
HeH'^-l-D  channel  near  threshold,  where  impulsive  dynamics 
are  inhibited  due  to  inefficient  transhitional-to-interiKil  en¬ 
ergy  transfer.  At  high  total  energy,  the  two  sources  for  oppo¬ 
site  isotope  effect  balance  and  branching  ratios  are  close 
to  1. 

TWQS-CC  integral  cross  sections  predict  a  significant 
alignment  effect,  favoring  k^  reactants  for  i;=0-3  in  the 
HeH^-hD  channel  and  for  i;=0  and  1  in  the  HeD‘^-»-H  chan¬ 
nel.  The  alignment  effects  are  consistent  with  cla.ssical  steric 
arguments  favoring  a  collinear  approach  geometry,  angular 
momentum  conservation  constraints,  and  energy  transfer  dif¬ 
ferences  in  impulsive  dynamics.  The  calculated  7  and  energy 
dependent  reaction  probabilities  are  consistent  with  this  pic¬ 
ture  apart  from  significant  quantum  scattering  resonances. 
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